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Abstract. We study a general class of models whose classical Hamiltonians are given 
hy H = U {x)p + V{x) /p, where x and p are the position and momentum of a particle 
moving in one dimension, and U and V are positive functions. This class includes the 
Hamiltonians Hi = x(p+ 1/p) and Hu = {x + l/x){p+ 1/p), which have been recently 
discussed in connection with the non trivial zeros of the Riemann zcta function. We 
show that all these models are covariant under general coordinate transformations. 
This remarkable property becomes explicit in the Lagrangian formulation which 
describes a relativistic particle moving in a 1+1 dimensional spacetime whose metric 
is constructed from the functions U and V. General covariance is maintained by 
quantization and we find that the spectra are closely related to the geometry of 
the associated spacetimes. In particular, the Hamiltonian Hi corresponds to a flat 
spacetime, whereas its spectrum approaches the Riemann zeros in average. The latter 
property also holds for the model Hu, whose underlying spacetime is asymptotically 
flat. These results suggest the existence of a Hamiltonian whose underlying spacetime 
encodes the prime numbers, and whose spectrum provides the Riemann zeros. 



PACS numbers: 03.65.-w, 02.30.Tb, 03.65.Ge, 03.65.Sq 
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1. Introduction 

In 1999 Berry and Keating conjectured that an appropriate quantization of the classical 
Hamiltonian H = xp, of a particle moving on the real line, could provide the long sought 
spectral realization of the Riemann zeros PQ [2] . These authors were led to this idea by 
the similarity between the semiclassical spectrum of a regularized version of the xp 
model and the average distribution of the Riemann zeros. The regularization introduces 
the constraints > ix and \p\ > ip in position and momentum, such that the product 
of their minimal values is equal to the Planck constant [ixip = 27rh). This proposal 
was made in the framework of Quantum Chaos and spectral statistics [21 HI [3] . About 
the same time, Connes proposed another regularization of xp based on the constraints 
|x| < A and \p\ < A, where A is a common cutoff [H]. In the limit where A is sent 
to infinity one obtains a continuum spectrum where the Riemann zeros are absortion 
spectral lines, according to Connes. This interpretation underlies the adelic approach to 
the Riemann hypothesis. These results have motivated several works in the last years 
on the xp model, and related quantum mechanical models, for their possible connection 
with the Riemann zeros [Tj-[2I] (see [22] for a review on physical approaches to the 
Riemann hypothesis). 

Specially relevant to this paper are the recent works [HI [19], which propose 
two different modifications of the xp Hamiltonian in order to have bounded classical 
trajectories and a discrete quantum spectrum. In reference [18], the classical 
Hamiltonian is Hi = x{p + ip/p), which adds to xp a non standard term xi^/p, where 
ip is a constant. The latter term implements, in a dynamical way, the constraint 
\p\ > ^p, but one still needs the constraint x > ix- The classical Hamiltonian Hi can 
be quantized in terms of a self-adjoint operator whose spectrum agrees asymptotically 
with the first two terms of the Riemann-Mangoldt formula that counts the number of 
Riemann zeros [23j . The Hamiltonian Hi breaks the symmetry between x and p, which 
is an appealing feature of the xp model. This fact led Berry and Keating to propose a 
new Hamiltonian Hu = {x + il/x){p + ip/p), which restores the x — p symmetry and 
implements dynamically both constraints on x and p, as can be seen from the appearance 
of the constants ix^p in it. 

The aim of this paper is to generalize the previous models, considering Hamiltonians 
of the form H = U{x)p + V{x)/p, where U{x) and V{x) are positive functions 
defined on intervals of the real line. This class of Hamiltonians have the remarkable 
property of being general covariant, which means that they maintain their form 
under general coordinate transformations, i.e. diffeomorphisms x' = f{x). These 
transformations change the functions U and V, according to prescribed laws, but 
the physical observables, such as energies, remain unchanged. General covariance is 
a signature of gauge symmetry, as it occurs in General Relativity. Indeed, we shall 
show that the present models describe the motion of a relativistic particle moving in a 
1+1 dimensional spacetime whose a metric can be constructed in terms of the functions 
U and V. The classical trajectories of the Hamiltonian H = U{x)p + V{x)/p, being 
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the geodesies of that metric. Hence these generahzed xp models acquire a geometrical 
interpretation which gives new insights into their quantum properties, and in particular 
their spectrum. 

The organization of the paper is as follows. In section |2] we introduce the classical 
models and show their general covariance. In section |3] we pass from the Hamiltonian 
to the Lagrangian formulation and present a relativistic spacetime interpretation, which 
is illustrated with several examples. In section |4] we discuss the classical trajectories in 
the Hamiltonian and Lagrangian formulations. In section |5] we analyze the semiclassical 
spectrum of the models introduced in section [3j We quantize the models in section [6] and 
show that general covariance is maintained. Finally, we present our conclusions. We have 



included in [Appendix A| the derivation of the inverse of the semiclassical quantization 



formula, and in 



Appendix B the quantization of the Hamiltonian H = p + li/p 



2. The classical Hamiltonian 

Let us consider a general class of Hamiltonians of the form 

Vix) 

H = U{x)p + xeD, (1) 

where x and p are the position and momentum of a particle moving in an interval D 
of the real line, and U{x) and V{x) are positive functions in D. We shall be mainly 
concerned with intervals that are halflines, D = {i^, oo), and eventually with segments 
i.e. D = {£x, £x), The two examples discussed in the introduction correspond to [HI [19] 



Hi =x(^p+^-^y D = {lx,oo) (2) 

Hn= (x + ^j^(^p+^Y D = (0,oo). (3) 



Berry and Keating also studied the model (|3j) on the whole real line, but we shall not 
consider this case here because the corresponding functions U and V are not positive. 
The positivity conditions on U and V are necessary, in order to have bounded classical 
trajectories, but not sufficient, as shown by the example H = p + i^/p (see section |3] 



and Appendix B ) . It is convenient to write U and V as 

U{x) = u'^{x), V{x)=v'^{x), (4) 

where u{x) and v{x) will also be positive functions. The Hamiltonians ([T]) change their 
sign under the time reversal transformation, i.e. 

X — )■ X, p — 7- —p =^ H — 7- —H, (5) 

which implies that if {x{t),p{t)} is a classical trajectory with energy E, so is 
{x{t), —p{t)} with energy —E. Upon quantization, the spectrum will contain time 
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conjugate pairs —En}, for appropiate boundary conditions related to the self- 
adjoint extensions of ([T|. The breaking of the time reversal symmetry is suggested 
by the statistical properties of the Riemann zeros, that are described by the Gaussian 
Unitary Ensemble distribution (GUE) [25]. 

The Hamiltonian ([T]) is covariant under general coordinate transformations of the 
variable x. Indeed, let us consider the infinitesimal canonical transformation 

x' = x + e{x), p' = {1 — dxe{x)) |e(x)| << 1, (6) 

that preserves the Poisson bracket 

{x,p} = l^{a;',j)'} = l + 0(e^). (7) 

Substituting these eqs. into Q, one obtains 

H{x, p) = [U{x') - e{x')d^U{x') + {d,,e{x'))U{x')] p' (8) 

+ \V{x') - e{x')d,V{x') - {d,,e{x'))V{x')] + 0{e^) 

p' 

= H{x',p') + Oie'), 
which has the same form as ([T]), for redefined functions 

U'ix') = U{x') - t{x')d.JJ{x') + (9,,e(x'))f/(x ), (9) 
V'{x') = V{x') - e{x')d,Vix') - id^,eix'))V{x'). 

These equations are the infinitesimal version of the transformation laws of one 
dimensional tangent and cotangent vectors 

U'ix') = ' Uixix')), V'ix') = %V{x{x')). (10) 

The momentum p transforms also as a cotangent vector (i.e. one form). Another way to 
state these transformations laws is by saying that the products U {x){dx)~^ , V{x)dx,pdx 
and also u{x){dx)^^^'^,v{x){dxy^'^ are invariant under reparametrizations of x. To 
preserve the positivity of the new functions U' and V, we shall restrict ourselves to 
diffeomorphisms x' = f{x), such that df{x)/dx > 0. These diffeormorphisms form the 
group denoted as Diff^. The interval D is mapped into the new interval D' = f{D). 
All the models related by diffeomorphisms are equivalent at the classical level. We shall 
organize them into equivalent classes described by the quotient 

Mci = {U,V,D}/I)m-^. (11) 

Each class can be uniquely characterized by a Hamiltonian which has a particularly 
simple form, 

w{x) = U{x) = V{x) ^ H = w{x) (^p + ^ . (12) 
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Any other Hamiltonian can be brought into this form by a convenient reparametrization. 
For example, the models ([2]) and ([3]) correspond to 



Hi — )■ wi{x) 
Hii wiiix) 



X 



X 



D = {h,oo) 
D = (0,00 



h = ejp (13) 
h = ejp (14) 
X in both cases. The function w(x) 



where we made the change of variables x' = 
is unique, up to the shift x — )■ x + cte. We shall call the canonical form (12) the 
symmetric gauge. Other gauges are possible, as for example U{x) = 1, which will be 
briefly discussed at the end of the next section. w{x) is a scalar function that can be 
computed in any coordinate system as 

w{x) = u{x) v{x). (15) 

Using the transformation laws of u{x) and v{x) one can verify that w'{x') = w{x), as 
claimed above. To find the coordinate transformation that brings a model into the 
symmetric gauge consider the equation 



v'{x') 



dx' 



v{x) 



u'{x') u{x 
In the symmetric gauge u'{x') ■ 
i.e. 



dx. 



(16) 



v'(x') 



so integrating (16) yields the mapping x' = f{x) 



x' 



4 = /(^) 



dy 



v{y) 



which is invertible, x = f 



ix' 



uiy) 
since v{x)/u{x) 



> 0. The constant 



(17) 
is left 



undetermined by this map, so it can be choosen at will. w'{x') is obtained using eq.(15) 



and the inverse of (17) as 



w'{x') = w{x) = u{x) v(x) = u{f '^(x'))v{f ^{x')). 



(18) 



An application of equations (17) and (18), is to show that apparently different models 



may turn out to be equivalent, as shown by the following case. Consider the model. 



,2P 



X 



(0,oo) 



Hiu = xp + ii- + it-, D 
X ^ p 

which, as the model II, is symmetric under the interchange x 



(19) 

p, differing from it in 



the term i^ip/xp in the Hamiltonian. The map (17) becomes in this case 

" Jo ^ Vv^x 



X 



X 



which plugged into (18) yields 



w 



III 



{x') = ^pv^^^T^ = x'. 



(20) 



(21) 



so that this model actually coincides with the model I defined in (13). 



In the definition of the family of Hamiltonians ([T]) we have imposed the positivity 
condition on U{x) and V{x). Let us suppose for a while that V{x) = 0. One can see 
that a reparametrization x ^ x' can bring U{x) to x' and so, all the Hamiltonians of 
the form U{x)p, are equivalent to xp. 
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3. Lagrangian formulation: relativistic spacetime picture 

An essential feature of General Relativity is that the fundamental equations of the 
theory take the same form in all coordinate systems. As shown in the previous section, 
this is also a feature of the models defined by the Hamiltonians ([T|, with respect to the 
coordinate x. Henceforth, one may suspect the existence of a general relativistic theory 
lying behind the models ([T]), which would provide them with a spacetime interpretation. 
In this section we shall show that this is indeed the case via the Lagrangian formulation. 
The Lagrangian associated to the Hamiltonian ([T]) is given by 

Vix) 

L = px- H = px-U(x)p —. (22) 

p 

In standard classical mechanics, the Lagrangian can be expressed solely in terms of x 
and X = dx/dt, a.s L = mx^/2 — V{x), where m is the mass of the particle and V{x) is 
the potential. To find L{x,x) in our case, we use the Hamilton equation of motion 

OH , V(x) , , 

to eliminate p in terms of x and x. This gives two solutions 

that depend on the sign of the momenta, r] = signp, which is a conserved quantity. The 
positivity of U{x) and V{x), imply that the velocity x must never exceed the value of 




U{x), for the momentum not to become an imaginary number. Substituting (24) into 



(22), yields a Lagrangian, 



L^(x,x) = -2ri^/V{x){U{x) - x). (25) 



for each value of rj. Notice that eq.(24) is singular if V{x) = 0, so that the Lagrangian 
cannot be expressed in terms of x and x. This is precisely the situation of the usual xp 
Hamiltonian, whose Lagrangian, L = p{x — x), has to be considered as a function of the 
three variables x, x and p. Later on we shall give an interpretation of this peculiar fact. 

At the classical level we can restrict the motion of the particle to a definite value 
of rj, but not at the quantum level, where both signs would be required. The action S 



corresponding to (25) is ( we choose rj = 1 



S = j dtLi = -'^ j \/U{x)V{x) {dtf - V{x)dtdx, (26) 

and it coincides with the action of a particle moving in 1+1 dimensional spacetime with 
metric g^^y, i.e 



u 



S = - I = - I ^/~g^ludxl^dx•', (27) 
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where a parametrizes the worldhne( we have set the mass of the particle to 1). Making 
the identifications 



X 







X 



X, 



the metric tensor becomes 

(700 = -AU{x)V{x) = -4W{x), goi = gw = 2V{x), 
In our conventions, the square of a hne element dx^ will be defined as 

{dsY = g^ydx^dx"", 



(2J 



^711 = 0. (29) 



(30) 



so that a time-like distance corresponds to {ds)"^ < 0, and a space-like distance to 
{dsY > 0. Eq.(29) imply that, under general transformations of the coordinate x, the 



function W{x) is a scalar, while the V{x) is a cotangent vector, in agreement with the 

(31) 



results of section [2j The determinant of the metric (|29|), given by 
g = det gf,^ 



-AV^(x) 



implies that g^j^^, is a non degenerate Minkowski metric since V{x) > 0. 

To gain further insight into the spacetime structure underlying the xp models, 
we shall employ the light-cone formalism, which we pass now to describe. Any two 
dimensional metric is conformally equivalent to a fiat metric. This means that it can 
be written as 



(ds)'^ = e^^ dx~^dx 



X = X[x ,x ), 



(32) 



where x^ are the light-cone variables and e-^ is the conformal factor. To find the 
transformation from the variables to the light-cone variables x~^,x~, we use the 

transformation law of the metric tensor 



9a 



dx^ dx'' 



\x 



X] 



dx"^ dx'^ 

where the ligh-cone metric corresponds to 



9++ = 9- 



0, 



9+- 



9- 



,2X 



(33) 



(34) 



Eqs. (29) and (34), allow us to write (33) as 



= [Vr(x^)9+x° - V{x^)d+x'^] (35) 

= \W{x^)d_x^ - V{x^)d^x^] (36) 

= - 8Vr(x^)5+x°9_x° + W{x^) + d+x^d^x^] . (37) 

Let us suppose, for a while, that x° depends non trivially on and x~, i.e. d±x^ ^ 0. 



Hence eqs. (35) and (36), would imply 

d±x^ 



f/(xi 



d. 



dy 

W) 



X 



dy 



+ cte. 



(38) 



General covariant xp models and the Riemann zeros 



8 



so that would be a function of x^, which is a contradiction since they are independent 
variables. We shall make the choice that x'^ only depends on x'^ 



X 



f{x' 



Eq.(36) is fulfilled automatically, and eq.(35) becomes 

d^x^ 



d^x' 



U{x^ 



dy 
U{y) 



X 



dy 
U{y) 



g{x 



and so 



dy 

U{y) 



f{x+)+g{x'), 



(39) 



(40) 



(41) 



where f{x~^) and g{x ) are generic functions. x^{x^,x ) is given, in an implicit way, 



by eq.(|41j). Finally, eqs.(37) and (41) provides the conformal factor, 

e^x = AW{x^)d+f{x^)d.g{x-). 



(42) 



Equations (39) and (41) give the map from x ' to x . However the map is not unique 
due to the freedom in choosing /(x+) and g{x'). This simply reflects the invariance of 
the metric (32) under general conformal transformations, x^ — )■ f{x^) and x~ — )• g{x~). 
In the conformal gauge (34), the tensors and connections simplify considerably. 

(43) 



The Christoffel symbols, F^^, have only ± non vanishing components 



vx^ = g^-d+g+^ = 2d+x, r:_ 

so that the equations of the geodesies a;^(s) read 



9 d-g+- = 29_x, 



2^± 



d^x 



ds^' 



+ 2 



dx 



ds 



d±x = 0, 



(44) 



where s is the propertime, i.e. {ds^ = 
= = -29+9_x, 
and the Ricci scalar, R = g^^R^j.^, 
R = -8e-2^(9+(9_x. 
It is not difficult to show that 



-e^^dx^dx . The Ricci tensor, -R^j,, becomes 



R 



++ 



R. 



R{x) 



Vix 



d^Wjx) 
V{x) 



(45) 



(46) 



(47) 



In the symmetric gauge (12), one has V{x) = w{x) and W{x) = w (x), so (47) becomes 



R{x) = -2-^>(") 



w{x) 



(48) 



We shall use this formula to relate different xp models to the underlying spacetime 
geometries. 
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In flat spacetimes the scalar curvature vanishes. Equation (48 ) provides the function 
w{x) corresponding to these cases 



R{x) 



0, MxeD 



w[x) 



ax + c, 



a > 0. 



(49) 



The condition a > comes from the positivity of w{x). If a > 0, the shift x ^ x — c/ a, 
brings w{x) to the form 



w{x) = ax, 



a > 0, D = {h, oo) 



(50) 



For a = 1, this model coincides with (13). The value of Wq = w{h) = ah, is independent 
of reparametrizations. If a = 0, w{x) is constant and the Hamiltonian is simply 



w{x) = c > 0, 



H = c\ p + 



P 



D = (0,oo) 



(51) 



where we have choosen the origin as the boundary of the interval D (see Appendix B ) 
Berry-Keating model 



This model was defined in eq.(14). The spacetime has a scalar curvature 

h^ 



w{x) = X + 



X 



R{x) 



Ah' 



(52) 



which is always negative, vanishes asymptotically as x ^, and diverges at the origin as 

2 



X 



Spacetimes with constant negative curvature 



Eq.(48) admits a solution with constant negative curvature 

R{x) = -\R\, \JxeD ^ w{x) = wocosh{x^/\R\/2), wq > (53) 
where > to guarantee the positivity of w{x). We shall take D = (0, oo). There 



is also a solution of (48) with positive curvature involving the cosine function, but it 



requires finite D domains in order to maintain the positivity of w{x). We shall not 
consider this case below. The interest of solution ( [53| is that the semiclassical spectrum 
coincides with that of the harmonic oscillator (see section [s]). 

Linear-log model 



An interesting variation of the linear potential (49) is to add a sub leading 
logarithmic term, i.e. 



w{x) = ax + /3logx =^ R{x) 



2/3 



x^(ax + /3 \ogx] 



2/3 



ax-" 



(54) 
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where the curvature decays asymptotically as x ^, with a sign determined by that of /3 
(a > 0). 

Power like models 

These models are defined by 

wix) = Ax" (A a > 0) Rix) = ^— — -. (55) 

We impose the condition a > to have a monotonic increasing function w{x). If a 7^ 1, 
the curvature vanishes asymptotically as x~^, and its sign is negative for a > 1 and 
positive for < a < 1. In section |5] we shall show that the asymptotic behavior of the 
curvature is intimately related to the semiclassical spectrum of the model. 



3.1. Flat spacetimes 



Let us study in more detail the model (50). Since the curvature vanishes, there is a 
choice of /(x"*") and g{x~) for which the conformal factor is constant, and therefore 
the geodesies are straight lines, it is given by 

1 



2a 



log 2;, 



which plugged into eqs.(39) and (41) yields {h = i^) 
1 



X 



2a 



logx 



X 



h{x X 



-^l/2 



and conformal factor (recall (42)) 
= h. 

The line element 

{dsY = h'^dx~^dx~, 



(56) 



(57) 



(58) 



(59) 



implies that the geodesies are straight lines in the x^ plane. Not the whole x^ plane 
is available for the motion of the particle because it is constrainted to the interval 
D = {h, 00). In light-cone coordinates the spacetime domain, U, can be obtained from 
eq.([57l) 



U = {(x+,x") |x^ e (0,00) 



x^x 



>!}■ 



(60) 



If x~^ and x~ denote the vertical and horizontal axes of the plane, then U is the 
region in the first quadrant that is above the hyperbola x"'"x~ = 1. This hyperbola is 
the worldline of the point x^ = h. More generally, the worldlines of any point x^ > h, 
are given by the hyperbolas x^x~ = [x^/h)"^. x~^ is the light-cone time coordinate and 
it fiows upwards. Eliminating x"*" in eqs.(57) one finds 



X 



he" 



X 



(61) 
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Figure 1. The region in shadow represents the universe U in light-cone variables 
defined in eq.(60|. Left: the hyperbola (dotted line) corresponds to the worldline of 
a given position = V x+x^ 

Right: 



cte, and the vertical (dashed) line corresponds to a 
light ray emanating at the boundary (eqs. 



(61)). 



the worldline of a particle 



with constants energy E, which bounces off at the boundary. 



Hence, the vertical lines, i.e. constant values of x , coincide with the classical solutions 
of the Hamiltonian, axp, namely x oc e°*. The line element (59) vanishes along these 



trajectories, which therefore represent light rays that start at a point on the boundary 
x^ = h and scape to infinity as x° — )■ oo (see fig. [T]): 

x~ = cte O light ray. (62) 



The line element (59) also vanishes along the horizontal lines, x~^ = cte, but they do 
not correspond to light rays since the time coordinate x^ is frozen. In this theory, the 
light rays are right movers. The left moving light rays are absent. This chirality is a 
reflection of the time reversal symmetry breaking of the Hamiltonian ([T]). 

The causal cone, i.e. {dsY < 0, at each point of U, is given by the second and 
fourth quadrants, which correspond respectively to the future and past events relative 
to that point. A particle follows straight lines, with negative slope that start and end 
at the boundary x^ = h. To show this fact explicitely, we solve the classical equations 
of motion of the Hamiltonian with w{x) = ax, for positive energy E (see section |4]) 

^2 ^ E^2a(t-to) _ ^4ait-to) ^ p2 ^ ^^-2a(t-to) _ ^ > g. (63) 

a a 



In light-cone variables (57) this equation becomes a straight line 



E 

g-" x+ + g° = — , E>0, q = e^^o/i^" > o, wq = ha, (64) 

Wo 

where q parametrizes the slope that depends on the time to where x(to) = b(^o)|- This 
line ends and starts at the points A and B of the boundary x~^x~ = 1, with coordinates 



xX^ = ^ = £- ( E ± ^E^ - Awl ) = e±-, (65) 
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where 



cosh(ae) 



E 

2wo 



> 1. 



(66) 



The energy E of the classical orbits are bounded by 2wq. The entire worldline of a 



particle with energy E, is given by a polygonal line made of linear segments (64), that 



come from the horizontal axis, x° — ?■ — oo, and approaches the vertical axis, x° — )■ +oo 
(see fig. [T]). The value of q, that parametrizes each segment, can be found matching the 
initial and final positions of consecutive segments, i.e. 



X 



Br. 



n 



— oo, oo 



(67) 



which means that the particle in the {n — segment bounces off sX x^ = h and starts 
a new orbit corresponding to the n^^ segment. This polygonal worldline represents a 
periodic motion, since after a shift x° — )■ x*^ + T^, i.e. 



X 



the equation 

(64) remains invariant if g — )■ e^"^^ q which, according to (67), gives the period as a 



function of the energy, i.e 

1 



cosh 



a 



E 

2wo' 



(6J 



Let us next study the model defined in equation (51), which also has a vanishing 



curvature. A choice that leads to a constant conformal factor is 

f{z) = g{z) = z, 



which using eqs.(39) and (41) yields 
„o 



X 



X 



X 



c{x + X ), 



and 



(69) 



(70) 



(71) 



= 2c. 

The constraint > provides the domain of spacetime 

U = {{x+,x-)\x+ + x- >0}. (72) 

which is depicted in figj2| which also shows the light-rays and the worldline of the points. 
The classical equations of motion have the solutions (we choose E > 0) 



x = c(t-to)(l-e^), p = e^ (73) 
x = c{t-to){l-e-'), p = e-^ (74) 



where e is defined as 



coshe = — , e >0. (75) 
2c 
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Figure 2. The region in shadow represents the universe U in Ughtcone variables 
defined in cq.(72). Left: the hyperbola (dotted line) corresponds to the world line of a 
given position = c{x'^ + x~) = cte, and the vertical (dashed) line corresponds to a 



light ray emanating at the boundary x'^ 



— 0. Right: worldline of a particle with 



constants energy describe by eqs.(76) and (77 1 with to, = 0. 



The solution (73) describes the particle moving towards the origin, i.e. x < 0, which it 
is reached at t = to. At that moment it bounces off and starts to move to the right, i.e. 



X > 0, as described by eq.(74). In the lightcone coordinates (70), eqs.(73, 74) becomes 



+ e'/^x- 
+ e'/^x- 



2to sinh(£/2), 
2to sinh(£/2). 



(76) 
(77) 



Fig. |2] illustrates the form of these trajectories. 

A conclusion of the results we obtained in this subsection is that the scalar curvature 
by itself does not fully characterize a model and that the boundary of spacetime may 
play an essential role. 

Before we leave this section we want to make some remarks: 

• The reformulation of the Hamiltonians ([T| as relativistic models described by the 



Lagrangians (25) holds strictely speaking at the classical level. At the quantum 



level, this relation is more involved. Indeed, one should start from the path 
integral in phase space and perform the integration over the variable p. For 
Hamiltonians of the form H = p^ /2m + V{x), this integral is gaussian and one 
gets the familiar Feymann path integral with Lagrangian L = mx^/2 — V{x). 
However for the Hamiltonians ([T]), the integration over p is not gaussian and in 



the saddle approximation, one gets extra terms in addition to the Lagrangian (25). 



In any case, the quantization of these models will be done in section [6] using the 
Hamiltonian formulation, where this issue does not arise. 



Given a model in the symmetric gauge (12), one can find a new coordinate x' such 
that U{x') = 1, that we shall call the p— gauge. The transformation x' = f{x), and 
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the new function V'{x') = l^(x'), can be derived from eq.(lO) and the scalar nature 
of W{x) 

dy 



x'-i', = f{x) 



w{yy 



v,{x') = w\r\x)). 



In the p— gauge the Hamiltonian takes the canonical form 



H = p + 



P 



whose square 



H'=p' + 2V'{x) + 



(7J 



(79) 



(80) 



coincides with the standard Hamiltonian with a positive potential 2Vp{x), except 
for the Vp{x)/p'^ term. One finds for example that the functions Vp{x) associated 



to the models (13) and (14) are given by 



2 2x 



Hi Vni(x) = h'e 



2„2x 



Hii V^iiix) 



D = (0,cx)) 
D = (0,oo) 



(81) 
(82) 



1 - ' 

which seem to have some relation with the Liouville model and the Morse potentials 
studied in [13] (see references therein), whose spectrum is related to the Riemann 
zeros in average. 



4. Classical trajectories and equations of motion 



In the symmetric gauge (12), the classical trajectories in phase space are curves with 
constant energy E, 



E = w{x) I p + 



P 



(83) 



For each position x, there are in general two different values of the momentum p that 
we denote as 



p±(x, E) = ^ [e ± ^E^-Aw^ix)) , 
and satisfy the relation 

1 

p+{x,E) 



P-{x,E)' 



^5) 



In fig. 3 we plot the classical trajectories corresponding to the models (13) and (14). In 



all these models, the trajectories are clockwise for positive energy and anticlockwise for 
negative energy. This is a consequence of the time reversal breaking of the Hamiltonian 
([T]). In the model (13), the particle starts at x = h, with a high momentum, say Phigh- 




Figure 3. Classical trajectories of the Hamiltonians ^ (left) and ([3| (right) in 
phase space with E > 0. We include for comparison the classical trajectory of the xp 
model, given by the hyperbola E = xp (dashed lines).. The parameters of the latter 
Hamiltonians are choosen as = £p — 1. The classical trajectories with negative 
energy can be obtained from the trajectories with positive energy replacing p — > —p. 



Then it moves to the right, while its momentum decreases until \p\ = 1, where x reaches 
a maximum xu{E). Afterwards, the particle moves back towards lower values of x 
and \p\. When the particle reaches the boundary of the interval, i.e. x = h, with a 
momenta |piow| < 1, it bounces off acquiring the original momenta phigh 



satisfy eq.(85). The latter relation can be written as log \p 



high I 



- log \P] 



low 



1/piow, that 
which is the 

analogue of the reflection of a particle in a wall. After the reflection, the particle restarts 
the motion from its original position and momentum, so that the classical trajectory is 
closed and periodic 



In the model (14), there is no a wall for the motion of the particle which follows 
a continuous and differential orbit around the point [h, 1] 



in phase space [19]. In both 
models, for large values of x and p, the classical trajectories approach the hyperbola 
E = xp, but they departure from it whenever x or \p\ are of the order of the parameters 
h or 1. 

The Hamilton equations of motion for the general Hamiltonian ([T]) are given by 



X 



wix] 



p 



—W [XI 



p + 



p 



(86) 



where x = dx/dt, w'{x) = dw{x)/dx. Computing x, and expressing p as a function of 
X and w(x), one finds 

6x 



X 



WiXjW [X) 



-4 + 



wix] 



w'^ix] 



(87) 



In section pi we formulated our model as a general relativistic theory. Equations (44) 



for the geodesies are therefore expected to follow from eq. (87). Indeed, let us choose 
for simplicity f{z) = g{z) = z, so that eqs.(39), (41) and (42) become 

w{x), 



X 



d±x 



.2X 



4w' 



X 
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where x^ = t and x^ = x. The proper time s and the time t along the trajectory are 
related by 

dt dx' 



so 



(dsy = -e^^dx+dx- = Aw\x) {dtf - 



w{x) 



dx 



X 



-1/2 



ds 2w{x) \ 
Taking a derivative respect to t in this equation yields 



(89) 
(90) 



d^ 
~d? 



dx~^ 
ds 



w[x 



w\x)x ^ w\x){xY ^ 



X 



w{x) 2w'^{x) 2w{x) 



which together with (87), leads to 

2 



2^+ 



d^x 



di' 



+ 



dx'^ 
ds 



2w'{x) = 0. 



This equation coincides with (44) (for follows from 

djf-w dxW 



d+X 



w'{x). 



w w 

The geodesic equation for x~ can be derived in a similar way. 



(92) 



(93) 



In the model (50) the solution of the eqs. of motion (86) is given by 

E 



^2 ^E^2a(t-to) _^4ait-to)^ p2 ^ ^^-2a(t-to) _ ^ > g, (94) 

a a 
where to is the instant where x = p. At the initial and final times, tij the particle is at 
X = h, which implies 



,2«(tj,,-io) 




SO that the period of the trajectory coincides with (68) 

E 



tf - ti 



— cosh ^ 

a 2wq 



1 E 
-log — , 

a Wo 



{E » Wo). 



(95) 



(96) 



The trajectories of the Berry-Keatig model with Hamiltonian Hu, (recall eq. (14)), 
are given in terms of elliptic functions [19]. In fig |4| we plot x{t) and p{t) for the 
Hamiltonians Hi and Hu with the same energy. The discontinuity of the curve for the 
Hi model is in contrast with its continuity for the Hu model. Despite of this fact, both 
curves have almost the same period. It is remarkable how fast the particle retraces its 
trajectory near the origin, which is the reason why the periods Tg in the two models 
converge for large trajectories. 

The general expression of the period of a trajectory with energy E, can be obtained 



integrating eqs. (86) 
Te = 



dx 



E 



E = 2w{xm) = 2w{xm)-, 



(97) 



w{x) ^E^-4w^{x)' 

where x^ and xm are the turning points. The minimal value Xm can actually coincide 
with the boundary value of the interval D, as for the Hi model. 
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Figure 4. Plot of x{t) and p{t) for the Hamiltonians Hi (continuous red curve) and 
Hii (dotted blue curve). 



5. Semiclassical analysis 



The number of semiclassical states with energy between and > 0, denoted as n{E), 
is given by the area in phase space swept by the closed trajectory, measured in units of 
the Planck constant 27rh. In the symmetric gauge it reads 

dx 



n{E) 



2Txh 



dx dp 



0<H<E 



2'Kh 



W[X] 



(98) 



where Xm,M are the turning points of the trajectory (recall (97)). We have not included 



in (98) the Maslow phase. The counting formula for the negative energy states is 
also given by (98) with \E\ = 2w{xm,M)- In the next section we shall see that this 



symmery is broken in general by the quantum model. In other gauges, the corresponding 



formula for n{E) is obtained from (98), by the replacements dx/w{x) — dx/U{x) and 
w'^{x) — )■ U{x)V{x), which shows that n{E) is invariant under reparametrizations of x 
(see section [2]). The derivative of the area of the trajectory, 27ihn(E), respect to E, 



gives the period (97). In the cases where w{x) is an invertible function, and assuming 
that Xra = Xq, one can invert eq.(98) in order to find the function w{x), or rather x{w), 



that produces a given n{E). The formula is given by (see Appendix A for the proof) 



X[W) — Xq 

2nw 



dEE- 



d [n{2E)\ 



wo 



dE 



E J ^w^-E^' 



(99) 



Let us apply eqs. (98) and (99) to the models studied in section 3 



Linear potential 



In the case of the linear potential 

w{x) = ax, D = {h, 00) 

one finds 



(100) 
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n{E) 




(101) 



E»l E 



2'Kha 



(\og—-l]+0{E-^) {E»wo) 
\ Wo J 



where Wq = w{h) = ah. The leading term, 0{E\ogE), depends only on a, and the 



next to leading term, 0{E), depends on a and wo- Let us compare eq.(lOl), with the 
Riemann-Mangoldt formula up to a height t 



{nit))^^(\og^-l)+l + Oit-'] 

ZTT \ ZTT / 8 



t » 1. 



(102) 



The first two leading terms in this formula agree with those of (101) under the 
identifications 

E 



t 



ha ' 



h = 27ch. 



(103) 



Setting a to 1, we recover the model I. Notice that the constant term in (101) vanishes, 
unlike in Riemann's formula where it is given by 7/8. 

Linear -log model 

One may ask which modification of w{x) would yield a counting function n{E) 



similar to (101), but containing a constant term, i.e. 

E 



n{E) 



27ih 



cosh ^ 



E 
2wo 




(104) 



To find this potential we use eq. ( 99 ) 



X — Xq = w — Wo + fih log 



1 - y/l - jWo/wy 

1 + - [wo/wy 



^ W — Wq + 2fih log 

whose inverse is 

w{x) ~ X + 2/i/ilogx, 



Wq 

2w'' 



X » 1. 



W » Wq, 



(105) 



(106) 



Based on this result we expect that a function w{x) reproducing the Riemann zeros 
would contain a log x in its asymptotic expansion. 



The Berry-Keating model 



The semiclassical spectrum associated to the function w{x) = x + h'^/x, defined in 
the halfline D = (0, oo), is given by [IH] 
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n{E) 



E»l 



E 






2'Kh 


K 


(■ 










E 


2-Kh 






1 



16/^2 \ / l6/i2 
2h^ log E 



(107) 



1 



where K and E are the elliptic integrals 

f7r/2 



T^h E 



+ .. 



Jo Jo 
Choosing t = E/h and h = 27rh, as in ( 103[ ), one gets an agreement with the first two 



terms of the Riemann-Mangoldt formula (102), in the large E limit. The constant term 



is absent in (107) which we believe is due to the absence of the term logx in w{x), 
for large values of x. All these features are shared with the linear potential analyzed 
previously. 

Spacetime with constant negative curvature 



The semiclassical spectrum associated to the function w{x) given by (53) and a 
domain D = {0, oo), is 

E 



E„ 



(109) 



where n and u are related to the parameters Wq and the scalar curvature R = — \R\ as 

It is an interesting fact that the harmonic oscillator spectrum is associated with a 
spacetime with constant negative curvature, at least semiclassically. This result does 
not contradict the chaotic spectrum associated to models defined in two dimensional 
euclidean spaces with constant negative curvature, since they have different spatial 
dimensions, i.v. one versus two 

Power like models 



To study the semiclassical spectrum of the models defined in (55), we shall 
distinguish the cases a < 1 and a > 1, which have very different properties. For 
functions w{x) that grow slower than x, one obtains 

< a <1 ^ n{E) (X E^, E^ oc n" (111) 

The density of states dn/dE behaves as Ea~^ and diverges in the limit a — > 0, where 
it becomes a continuum. The model a = corresponds to a constant w{x) and it has 



a continuum spectrum (see Appendix A). For potentials that grow faster than x, one 
finds 



a > I 



n{E) oc E + 0{Ei 



En oc n + 0{n") 



;ii2) 
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which is a hnear spectrum with corrections. In the hmit a — )■ oo, one recovers the 
harmonic oscillator spectrum, which is consistent with the exponential growth of w{x) 



(recall eq.(53)). 



In table 1, we collect the results obtained in this section and a comparison with 
the curvature of the associated spacetime models. From these examples, we draw the 
following conclusions: 

• In models where 'w{x) becomes linear for x >> 1, the first two leading terms of n{E) 
are of order ElogE and E. The converse is also true, i.e. the latter asymptotic 
behaviour of n{E) forces w{x) to be linear for x » 1. Correspondingly, the scalar 
curvature R{x) vanishes faster than 

• The appearance of a constant, in the next leading correction to n{E), requires a 
logaritmic term in w{x), in addition to the linear term. 

• Based on these arguments, one is lead to conjecture that the function w{x), whose 
quantum spectrum yields the exact Riemann zeros, is of the form 

w{x) = X + fi log X + ws_{x), (113) 

where Wfl(x) represents the fluctuation part of w{x). The role of Wfl(x) is to provide 
the fluctuation term nfi{t) of the Riemann-Mangoldt formula for the exact position 
of the Riemann zeros, 



nnit) = {n{t)) + ?2fl(t), ?2fl(t) = ^ argC Q + it 



where ({s) is the Riemann zeta function. It is not clear at the moment how to 
construct waix), or even if it exists. 



Model 


w{x) 


n{E) 


R{x) 


Linear potential 


X 







Berry-Keating potential 


X + h'^x'^ 


^\og^^+0{\ogE/E) 




Linear-log potential 


X + 2fi log x 


f log^+/. + 0(i?-i) 


4fix~^ 


Harmonic oscillator 


wq cosli(x/2/i) 


li{E/2wQ - 1) 


-1/2^2 


Sublinear potentials 


(0 < a < 1) 




-2a{a - l)x-2 > 


Superlinear potentials 


x" (a > 1) 


E + 0{E^I'^) 


-2a{a - l)x-2 < 



Table 1.- Semiclassical counting function n{E), and scalar curvature R{x), for the models 
discussed in the text. We show the asymptotic behaviour of the functions n{E) and R{x), 
i.e. E » l,x » 1, with h = l. 
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6. Quantization 

To quantize the classical Hamiltonian ([T| we shall choose the normal ordering 
prescription 

H = u(x) pu(x) + v(x) - v(x), X ED = (ix, oo) (114) 

P 

where p = —ihd/dx and 1/p is the one dimensional Green function 

{x\^\y) = -'j.0{y-x), (115) 

written in terms of the Heaviside step function. This normal ordering generalizes the one 
in references P^ [T^. One can check that 1/p is the inverse of p acting on wavefunctions 
ip{x) that vanish in the limit x — )■ oo, 

P T (a;) = - ^ ^ dy9{y - x)ij{y) = dy 5{y - x)ip{y) = %Ij{x), 
I \ f°° dip(y) 

- p i' ] (x) = — dy 9(y — x) — — = — lim ip(y) + i'M, (116) 

p J Ji,, oy y^°° 

where we have assume that x > l^- There exist other possible normal orderings defining 
H, as for example \ {v?{x)p+pv?{x)+v'^{x)p~^+p~^ v'^{x)) (see also pS]). However, the 



choice (114) yields a Schroedinger equation which, as we shall see below, is equivalent to 
a second order differential equation, suplemented with a non local boundary condition. 
Our construction parallels in that way the Schroedinger equation arising from the 



quantum Hamiltonian H = p^/2m + V{x). The action of (114) on a wavefunction 
il){x) is given by 

iH^Mx) = -ihu{x)^{u{xU{x)} (117) 

dx 

POO 

-ih~^ dyv{x) 9{y - x)v{y)ij{y). 



This operator will have a real spectrum if it is self-adjoint, which requires, first of all, 
that it be symmetric 



{tPi\Hij2) = {HHi^2), VV>i,V'2 e V{H), (118) 

where 'D{H) is the domain of the operator H. To study this condition, let us define the 
quantity [29] 



Using (117) one finds 



(119) 
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(120) 

ux 

c 

/'OO /»oo 
/ dx dy v{x)v{y)il)\{x)il)2{y){9{x — y) + 9{y — x)) 

/•OO poo 
/ dxdyv{x)viy)ri{x)ilj2iy), 

where we have assumed that hmj.^00 u{x) il'i,2{x) = 0. Hence ^ is a symmetric operator 



iff Q12 = which, in view of eq.(120), is guaranteed if iIji^2 satisfy the equation 



e'^hu{l^)%l){l^) + h-^ dxv{x)%lj{x) = (121) 
where d parameterizes the selfadjoint extensions of H. The Schroedinger equation of 



the Hamiltonian (117) is given by 



d 

-ihu{x)— {u{x)ip{x)}—ih^^ I dyv{x) 0{y—x)v {y)ip {y) = Eip{x), {122) 
dx IP 



which we write as 

u{x) d 



h 



iE 



{u{x)i){x)}-—-i){x) + h-^ / dy e{y-x)v{y)i){y) = Q. {12^) 



v{x) dx v[x) 
Taking a derivative respect to x and letting x = one gets 



d 

dx 

h 



u{x) d u \ 

a , , — {u[x)ip[x)\ ——ip[x) 



v{x) dx 
u{x) d 



v{x) 



h ^v{x)il){x) 



0,(124) 



v{x) dx 



{u{x)ip{x)} 



iE 



V[X] 



ip{x) 



+ h'' / dyviy)i;{y) = 0.(125) 
Jix 



It would seem that one has to solve simultaneously eqs.(124) and (125). However, for 



wavefunctions that decay sufficiently fast at infinity, the latter equation follows from 
the integration of the former one, dropping a term at infinity. Hence, the Schroedinger 



equation (122) is equivalent to the second order differential equation (124), which 



justifies the normal ordering (114). 



Eqs. (124) and (125) exhibit the general covariance discussed in the section 2 This 



can be shown using the transformation laws of u{x),v{x) given in section |2| together 
with that of ip{x), 



i:\x') {dx'f^ = tP{x) {dxf^^. 



(126) 



which preserves the form of the scalar product of the Hilbert space under 
diffeomorphism. 



dx' (^'(x'))*^'(x') 



dx ip*{x) ip{x). 



(127) 
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As was done in the classical and semiclassical analysis, it is convenient to work in the 



symmetric gauge. Let us transform eq.(124) into that gauge. First we write it as 



h 



u{x) d 
v{x) dx 



fC^^^ {u{x)i;{x)] 
v[x) dx 



iE 



u{x)v{x) 



u{x)ip{x) 



— u{x)ip{x) 



0, 



and make the transformation x — )■ x' (17) to the symmetric gauge (i.e dx' 
dx v{x)/u{x)), 



h- 



d 

dx' 



h 



d<p{x') 
dx' 



iE 



wix 



— 0[X 



0, 



(128) 



where w{x') is given by eq.(18), and the function is defined as 

(129) 

Mx). Now we rename as and 



0(x') = u{x{x'))il){x{x')). 
Notice that 0(x) is a scalar function, i.e. (f)'{x' 



write eqs. (128) and (121), for an eigenfunction (j)E{z), with energy E, as 



dz 



iE 



dz 



E{Zo 



w[z^ 



^e{z) = 0, Z>Zo 



dz 4>e{z) 



0. 



The norm of the wave function iIj{x) becomes in this gauge 

""^ dz 



dx ijj*{x) ^jj{x) 



w{z) 



0*(z)0(z). 



(130) 
(131) 

(132) 



Let us next discuss the behaviour of the eigenfunctions under the time reversal 
transformation. The classical model has the property that if {x{t),p{t)} is a trajectory 
with energy E, then {x{t), —p{t)} is a trajectory with energy —E. Under time reversal, 
the quantum Hamiltonian ( 114[ ) changes sign, which leads us to expect a relation between 
eigenfunctions with energies E and —E. Indeed, taking the complex conjugate of 
eqs. ( p0p3l| ), and comparing them with those for an eigenfunction 0_£;(z), one finds 



^e{z) oc ^-e{z) 



= or TT 



(133) 



Hence, if ^9 = or tt, the spectrum displays the time reversal symmetry E -H- —E (this 
fact was already observed in [IS])- The difference between these cases resides in the 



existence of a zero energy state. The solution of equation ( 130 ) for E = is given by 

<f)E=o{z) = Ae-'/''' + Be'^''' (134) 

The exponential growing term e^^'^ will typically yield unnormalizable wave functions, 
so we only consider the decaying term. In this case the nonlocal boundary condition 



(131) becomes 



h{e'^ + l)e 



-h/h 



TT 



(135) 
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so, only for this value of is 00 = e an eigenfunction of the Hamiltonian with a 
norm given by 

dz 



(^o|^o)= / (136) 
which will be finite for a large class of models that includes the ones studied in sections 



|3jand[5j In summary, when t9 = 0, the spectrum of the Hamiltonian (114) contains time 
reversed pairs {En, —En}, excluding the zero energy, while if -i? = vr, in addition to the 
time reversed pairs, there is a normalizable zero energy state. Depending on the form 
of 'w{z), the spectrum may, or may not, contain a continuum part. 

6.1. The model I 

In the following subsection we shall apply the previous formalism to the model with a 
linear potential w{x) [18] 



w{z) = z, D = {zo,oo). (137) 

The eigenfunctions of the Hamiltonian ( |114 ) are the solutions of the differential equation 
(130), which in this case read {h = 1) 



nz) - -0'(z) + 



iE 



0. 



:i38) 



The solutions of this equation are given essentially by the modified Bessel functions, but 
only the iC-Bessel function gives a normalizable solution 

1 iE 

2 ~ Y' 



A,,z'-''Kjz) 



V 



(139) 



where is the normalization constant. Using eq.(129) one obtains 

ipEiz) = z^ Kl_^E{z), z > zo, 
whose asymptotic behaviour is 



ipsiz) oc 



Z 2 



1 , iE 

Z 2^ 2 e 



Zq < z « E/2 
z » E/2. 



(140) 



(141) 



which shows that in the region zq < z << E/2, the function "j/^e is given approximately 
by the eigenfunction of the Hamiltonian H = {xp + px)/2 [H |9l HO]. Notice that 
xm = E/2 coincides with the maximal elongation of the classical particle. Beyond this 
value the wave function decays exponentially, as corresponds to the particle entering the 
classical forbidden region x > xm- The equation for the eigenenergies can be obtained 
plugging (139) into (131), and using the integral 



dzz'-'KJz) 



K^^i{zo) 



(142) 
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with the result [18] 

The asymptotic behaviour of the fC-Bessel function 



0. 



(143) 



Ki 

2 



-7rt/4 



ze 



it/2 



t » 1 



yields the asymptotic limit of (143) 

fE fE\ i3 

cos — log — 

V2 ^V^oe; 2 

so that the eigenenergies E > behave as 



E » 1, 



n{E) 



E 
2^ 



log 



E 



2^ 



e z 



(144) 



(145) 



(146) 



The two leading terms in this equation agree with the Riemann formula (102) and 



the semiclassical result (101), with the identification zq 
eq.(103) (notice that wo = zq = h). We must set ^ = 



= 2tt, which is the same as in 
to guarantee time reversed 
eigenenergies in analogy with the symmetry of the Riemann zeros on the real axis. 
The quantization of the model brings in a constant factor —1/2 in the counting formula 



(146), so that the factor 7/8 of Riemann's formula remains unexplained. In this respect, 
we recall the comment made in the previous section that adding a term logo; to w{x), 
can give rise to this constant term in the counting formula. 



7. Conclusions 



In this paper we have studied the properties of a family of one dimensional classical 
models, and their quantized version, that are extensions of the well known xp model. A 
fundamental property of these models is that they are covariant under general coordinate 
transformations, so that they can be organized into equivalent classes that describe the 
same Physics. This fact suggests some some sort of universality that could perhaps be 
explained using renormalization group arguments, as the ones employed in reference [S], 
where the operator 1/p was also considered. 

General covariance manifests itself most clearly in the Lagrangian formulation of a 
relativistic particle moving in a 1+1 dimensional spacetime. The geometrical properties 
of this spacetime turn out to be related to the spectral properties of the associated 
quantum model in an deep way. To wit, when the curvature of spacetime vanishes fast 
enough at infinity, one obtains a spectrum that coincides with the Riemann zeros on 
average. It is tempting to think that certain fluctuations of asymptotically flat metrics, 
may cause fluctuations in the spectrum so as to accurately reproduce the Riemann 
zeros. It is expected that these fluctuations are determined by the prime numbers, but 
the precise manner in which this might occur is unclear. This general class of models can 
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arise as effective descriptions of the dynamics of an electron moving in the plane under 
the action of a perpendicular uniform magnetic field and a electrostatic potential of the 
form V{x,y) = U{x)y + V{x)/y, where x and y are the two dimensional coordinates 
[15] . Further investigation will be needed to clarify these issues. 

We hope the results presented here shed new light on the spectral interpretation of 
the Riemann zeros, stimulating the research into such a interdisciplinary topic 
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Appendix A. Semiclassical Abel like inversion formula 

The aim of this appendix is to derive a formula that, under certain conditions, permits 
to construct the potential w{x) that reproduces a given semiclassical counting formula 
n{E). For Hamiltonians of the form H = /2m + V{x), the answer can be obtained 
using the Abel inversion method [2H]- We shall also use Abel's method for Hamiltonians 
of the form H = w{x){p + First of all, we shall assume that w{x) is a monotonic 

increasing function in the domain D = {xq,oo), so that it is invertible, x = x{w). In 
this case the semiclassical formula (98) becomes 

1 r^' dx I 

n{E) = — ^E^-Aw\x), E = 2w{xm), (A.l) 



which we write as 

2TThn{E) 
E " 

Making the change of variables (with E > Q) 




dv 

XQ "'ro 



Differentiating with respect to s yields 



ds dr _ g2 



(A.2) 



/ N 1 2 1 , , , 

r[x) = — -, s = — , ro = — >s, wo = w{xo), (A. 3) 

w[x) E wq 



and using the inverse x{w)^ we transform (A.2) into 

f{s) = '7ihsn{s)= / dx a/ r2(x) — = / dr — Vr"^ — s^. (A. 4) 



dfis)_ r\^dx s 
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Next, we use the chain of identities 



ds 



1 



df{s) 



, r° , dx s 

Qg I (2f — — -^^^^^^^^^^^^^^^^^^ 
V/S2 - 1/2 ds Jy J, dr ^(s2-y2)(^2_32) 



(A.6) 



"° , dx r , s IT 

dr .A, v^(s2-y2)(^2_52) 2 



that can be obtained from the Frobenius theorem and the integral 



dx 



X 



TT 



^(x2-a2)(62_a;2) 



< a < 6. 



(A.7) 



Replacing f{s) = nhsn^s) into (A.6), and undoing the change of variables (A. 3), yields 
finally 



X[W) — Xo 

2hw 



d 

dEE 

wo 



dE 



n{2E) 
E 



1 



E^ 



(A. 



It is not guaranteed a priori that x{w) is an invertible function. This fact may restrict the 
functions n{E) that can be obtained as semiclassical spectrum. Another issue concerning 



(A.8) is the following. Adding to n{E) a term linear in E^ produces the same function 



as can be easily seen from (A.8). We lost track of this linear term when taking the 



derivative respect to s in eq.(A.4). To recover this term, one has to plugg the function 



w{x), obtained from (A.8), back into (A.l), and read the linear term in E. 



It is worth to compare these semiclassical formulas with those associated to the 
standard Hamiltonian H = + V{x). The analogue of eq. (A.l), for an even potential 
V{x), is given by 

2 



n{E) 



I dx ^JE-V{x), \E\ = V{x, 

JXQ 



(A.9) 



where we have ignored the Maslow phase, and the analogue of (A.8) is given by 

x{V) = hrdE^4^^=^. (A.10) 

JVo 



dE y/V -E' 



The latter equation was used by Wu and Sprung to obtain a potential V{x) whose 
semiclassical spectrum coincides in average with the Riemann zeros ([31], for a review 
see [22]), 

2 



x{V) 



TT 



V ^ (2V 
log 



Tre^ 



V{x) oc 



X 



log X 



X, V » 1 



(A.ii; 



This result was used as a seed for a numerical reconstruction of a potential whose 
spectrum matches a large number of Riemann zeros lying at the bottom part of the 
critical line. Quite interestingly, that potential has a fractal structure whose dimension 
is nearby the value 1.5. A problem with this approach is that the Hamiltonian is time 
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reversal invariant, a fact which does not agree with the distribution of the Riemann 
zeros which follow the GUE statistic which is characteristic of time reversal breaking 
random Hamiltonians (see [22]) for a discussion on this issue). 

Independently of the previous works, Mussardo employed (A. 10) to find a potential 



whose spectrum behaves, in average, like the prime numbers [32j. The prime number 



theorem (PNT) states that the number of primes up to x behaves asymptotically as 
7r(a;) ~ x/logx, so that their density decreases as dn{x)/dx ~ 1/logx, as conjectured 
long ago by Gauss and Legendre [23]. The PNT implies that the n*'^-prime number pn 
grows roughly as Pn ~ nlogn. Mussardo noticed that this growth allows one to find 
a quantum mechanical model whose energies are the primes numbers. Choosing the 
leading term in the expansion of the Riemann formula for 7t{x), 

7r(x)^Lt(x) = r ^ (A. 12) 

he obtained 



x(V) ~ =^ V(x) ~ (xlogx?, (x, V » 1) (A.13) 

log 1/ 

This result was also used as a seed for finding a potential whose spectrum matches 
precisely the lowest prime numbers [22]. As in the case of the Riemann zeros, the 
prime potential has a fractal structure with dimension near 2. The difference in fractal 
dimensions, 1.5 versus 2, is consistent with the fact that the Riemann zeros are less 
random that the primes numbers. The former ones satisfy the GUE statistics and the 
latter follow an almost Poissonian statistics. It is rather remarkable the proximity of the 
prime number /Riemann zeros potentials to the harmonic oscillator potential. In table 
2 we summarize these semiclassical results together with other well known cases. 



Model 


V{x) 


x{V) 


En 


n{E) 


Free particle 


cte 


< x < oo 


continuum 




Harmonic oscillator 




W 


n 


E 


Potential well 


cte 


< |x| < 1 




y/E 


Riemann zeros 


^ log X ' 


y/VlogV 


2iTn 
log n 


^los-^ 

2-K 27re 


Prime numbers 


(xlogx)^ 


Vv 
logy 


nlogn 


E 
logE 



Table 2.- Semiclassical spectrum associated to the classical Hamiltonian 

H = p'^ + V{x). 



Appendix B. Quantization of H = p + ip/p 

In standard textbooks of Quantum Mechanics it is taught that the momentum operator 
p = —ihd/dx is self-adjoint acting in the Hilbert space of square integrable functions, 
L'^{D), in two cases: i) D = R is the real line, and ii) D = {a,b) is a finite interval 
of the real line [281 EO]- In case i), the operator p is essentially self-adjoint, and in 
case ii) p admits infinitely many self-adjoint extensions characterized by the boundary 
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condition = e*'^^/'(a), where d G [0, 27r). However, p is not self adjoint when D is the 
halfline, D = (0, oo), and therefore its spectrum is not real. A solution of this problem 
is suggested by the model we have discussed in this paper. Indeed, let us choose the 
simplest non vanishing potential w{x), namely a constant 



WiXj 



H = p + 



P. 
p 



(B.l) 



defined on the halfline D = (0, oo). This model is equivalent to (51), via a scale 
transformation which gives the relation c = ip. We shall show below that H is self 
adjoint and that its spectrum is a continuum and eventually a bound state. This model 
illustrates in a simple example the more complicated models considered in the main 
body of the paper. In spite of its simplicity, this model shares several features with the 
so called Kondo model in Condensed Matter Physics, which suggests that it may have 
other applications. 



To quantize (B.l), we follow the steps of section |6| The Schroedinger equation is 
given by 



ih 



dip{x) 



dx h 

Taking one derivative respect to x, gives 



dydiy — x)tp{y) = Eip{x). 



,d^ilj{x) 



dx"^ 

whose general solution is 



+ iEh 



dip{x) 



dx 



+ £2V;(x)=0, 



ijE{x) = A{E) e 



ikj^(E)x 



+ B{E) e 



where 



k±{E) = j^{E± sign{E)^E^-Ul 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



We are assuming in (B.5) that E is real, for other values we replace sign(£') by ±E/\E\. 



Using (B.5) one can compute the von Neumann deffect indices n±, which give the 

(B.6) 



number of linearly independent solutions of the equation 

n± = dim {iIj± \H iIj± = ±izilj±, Imz > 0} 



Choosing z = 2ipC (c > 0), we get two normalizable solutions of (B.5), namely 



ip+{x) = y4exp 
ip-{x) = Aexp 



h 



(-C + + 1 



(B.7) 



= 1, 



hence n+ = n_ = 1, which implies, by the von Neumann theorem, that the operator 
H is self-adjoint, with an infinitely many extensions parametrized by the group U{1) 
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[28| |29] . The spectrum of H is given by two intervals whose boundaries are ±2^p and 
±00. and a bound state with eigenvalue Eq 



spec H = CU {Eo} = (-00, -2ip) U (2£p, 00) U {^o}- 



(B.8) 



We denote by C the continuum part of the spectrum. It is convenient to parametrize 
the two branches of the continuum as follows 



E = 2ipri cosh u, 



rj = sign(i?) = ±1, 
in which case the momenta (|B.5|) become 

k-{E) 



u> 0, 



h 



h ' 



and the wave function (B.4) 



(B.9) 



(B.IO) 



(B.ll) 



The discrete eigenvalue of H appears for \Eq\ < 2£p with a normalizable eigenfunction 
corresponding to the momenta i.e. 

\Eo\ < 2lp =^ i,E,{^) = Ce-'=°^ fco = -iK = ^ (-^^0 + ^4£2 - El) .(B.12) 

To find the value Eq, we impose the non local boundary condition (121), which 
guarantees that H is an hermitean operator, 



e'Vi.o(0) + f 



(B.13) 



where is the parameter that characterizes the self-adjoint extension of H. In eq.(121 ), 
we made the shift -i? — )■ -i? + vr, so that the first term in (B.13) changed its sign. Plugging 

(B.14) 



( |R12| into ( |R13| one finds 

-i-d 



kn 



which gives 



En 



h 



2£psint9, 



71 71 

2' 2 



(B.15) 



The restriction on -(9 comes from the relation cos'd oc Refco > 0. If 7r/2 < ^9 < vr, the 
equation (B.13) is not satisfied and therefore there is no a bound state. We shall restrict 
below to the case \^\ < n/2. If -i? = 0, one gets Eq = 0, so that the spectrum is time 
reversal symmetric. In the limits — )■ ±7r/2, one has Eq = ±2ip, and ko = ^iip/h, so 
that the eigenfunction (B.12), becomes a plane wave. The constant C in eq.(B.12), is 
fixed by the normalization of the wave function 

dx\^Eo{x)\^ = l^C = 



2ip cos ^ 



h 



(B.I6) 
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The measure of the size of the bound state is given by the average of x, 

h 



dx x\il) Eo{x)Y 



2ip cos 



(B.17) 



and diverges in the hmit i} — t- ±7r/2. The operator H is self adjoint, then the spectral 
theorem implies that its eigenfunctions ipE form an orthornormal basis, namely 



(V^EolV'iJo) 



{lpE\lpE') 



dxi!Eo{x)i'Eo{x) = 1, 

) 

dx'4)*E^{x)'4)Eix) = 0, 



EeC, 



poo 

/ dx^|J*E{x)^pE'ix) = 5{E - E'), E,E'eC. 
Jo 



(B.18) 
(B.19) 
(B.20) 



(B.18) coincides with (B.16). Eq. (B.19), gives the relation between the coefficients 



A{E) and B{E) of the wave function (B.ll) 



A{E) k*-ik+{E) e'^-ir]e 



B{E) kl-ik.{E) e'^ -ir]e-'''' 



(B.21) 



where we have used eqs.(B.lO) and ( |B.14 ). Condition (B.20), together with (B.21), fix 



the form of these coefficients. Using (B.ll) one obtains 



7r5{k'^ - k+) + iP- 



k' - k. 



+ Be' 



+ A%Be' 



n6{k'_-k+)+iP— 



1 



Be Ae' 



7r5{k'_ -k^)+iP 
1 1 



7r5(A;; - k^) + iP- 



k'~k. 



kL - L 

where k± = k±{E), k'_^ = k±{E') and denotes the principal part of K To derive this 
equation we have used the improper integral 

1 



dxe'"'' = 7T6{k)+iP 



k'' 



which is the integral version of the distribution identity 
— = -mSik) + p\. 



Eq.(B.20) is satisfied provided 



(B.22) 



(B.23) 



Al AE'6ik+ - k'^) + Bl BE'5{k^ - k'_) 
+ A*E BE'S{k+ - k'J) + B*E AE'6{k. - k'^ 



n 



-5{E-E'), 



and 



{A*j„B%)M 



Ae' 

Be' 



(B.24) 



(B.25) 
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where M is the matrix 



M 



1 1 



h 



k.-k', k.-k' 



Using 



- A;;) = ^ , " biE - E'), 



6{k± - k') = 0, 



^±1 



eq (B.24) becomes 



\A, 



I5f 



and similarly 



Asp e-" + 15^1^ e' 



2 „M 



27r;i 



Finally, eq. (B.21) implies 



(B.26) 



(B.27) 



(B.28) 



(B.29) 



A{E) 
B{E) 



^y 87r/i(cosh u — rj sin t?) 

i?7 e~" 



87r/i(coshM — ?7 sint?) 



(B.30) 



These expressions satisfy eq.(B.25), which finally proves that the wave functions 



(B.ll), with coefficients given by (B.30), together with the normalizable state (B.12) 



form an orthonormal basis. The operator if, has the physical meaning of momentum 



rather than energy. Its spectrum (B.8) almost coincide with that of the momentum 
operator p defined on the entire line, except in an interval around the origin (— 2£p, 2£p), 
which is replaced by a bound state localized at the edge of the system. This result is 
reminiscent of the Kondo model where a bound state is formed between an impurity 
localized at the origin and the conduction electrons [33J. However in our model there 
is no spin and it does not describe a many body system, so this analogy is for the time 
being formal. 
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